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The relationship between variability, luminosity and redshift in the South Galactic 
Pole QSO sample is examined in an effort to disentangle the effects of luminosity and 
redshift in the amplitude of the optical variations. The anticorrelation between vari- 
ability and luminosity found by other authors is confirmed. Our analysis also supports 
claims that variability increases with redshift, most likely due to an anticorrelation 
between variability and wavelength. In particular, our parametric fits show that the 
QSO variability-wavelength relation is consistent with that observed in low-luminosity 
nearby active galactic nuclei. 

The results are used to constrain Poissonian-type models. We find that if QSO 
variability results from a random superposition of pulses, then the individual events 
must have B-band energies between ~ 10 50 and a few times 10 51 erg and time-scales 
of ~ 2 yr. Generalized Poissonian models in which the pulse energy and lifetime scale 
with luminosity are also discussed. 
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INTRODUCTION 



O yariability is a conspicuous property of Active Galactic Nu- 
■ clei (AGN) over a wide range of luminosities (from LINERs 

t I to QSOs) and redshifts (from essentially up to z > 4). 
Given that the emissivity of QSOs peaks in the optical-UV 
range, creating the so-called Big UV Bump, the study of the 
variability properties in this spectral range is central to the 
understanding of the energy-generation mechanism. 

Determining the way in which the amplitude of the vari- 
ations depends on the source luminosity and redshift would 
be an important step towards constraining AGN models. 
Several observational studies have addressed this issue by 
using long term monitoring campaigns on large samples of 
objects, mostly carried out with photographic plates (An- 
gione 1973, Uomoto, Wills & Wills 1976, Bonoli et al. 1979, 
Trevese et al. 1989, Netzer & Sheffer 1983, Pica & Smith 
1983, Cristiani, Vio & Andreani 1990, Giallongo, Trevese & 
Vagnetti 1991, Hook et al. 1994, Trevese et al. 1994, Cris- 
tiani et al. 1996, Di Clemente et al. 1996). Despite some 
early disagreement, it is now firmly established that there is 
an anticorrelation between variability and luminosity, such 
that the relative amplitude of the variations decreases as we 
move up the luminosity scale (e.g. Hook et al. 1994, Cristiani 
etal. 1996). 

Given the strong correlation between luminosity and 
redshift in most samples, the relationship between variabil- 



ity and redshift is not easily disentangled from that between 
variability and luminosity. Cristiani et al. (1996) found a sig- 
nificant positive correlation between variability and redshift 
based on the analysis of 486 objects merged from three dif- 
ferent samples, in agreement with a similar trend previously 
reported by Giallongo etal. (1991). The latter authors, how- 
ever, did not find any significant correlation between vari- 
ability and luminosity, whereas Cristiani et al. did. A strong 
anticorrelation between variability and luminosity was also 
found by Hook et al. (1994) in an analysis of nearly 300 
QSOs in the South Galactic Pole region (hereafter the SGP 
sample). While Hook et al. did not rule out the possibility of 
a positive correlation between variability and redshift in the 
SGP data, they concluded that the observed correlation was 
consistent with being entirely the result of the L-z correla- 
tion propagated through the variability-luminosity relation. 
Though on balance the evidence seems to point towards a 
positive variability-redshift correlation, there is as yet no 
consensus (see discussion in Hook etal. 1994). 

The relationship between variability and redshift is a 
key link between the variability properties of QSOs and of 
less luminous AGN. IUE observations of Seyfert 1 nuclei and 
low-redshift QSOs have demonstrated not only that vari- 
ability is anticorrelated with luminosity, but also that the 
variations get larger towards shorter wavelengths (Edelson, 
Krolik & Pike 1990, Kinney etal. 1991, Paltani & Cour- 
voisier 1994, Di Clemente etal. 1996). If this behaviour also 
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holds for high-redshift QSOs, one would expect variability 
studies carried out at fixed observed wavelengths (say, the 
B-band) to show a positive correlation between variability 
and redshift (Giallongo et al. 1991, Hook et al. 1994, Trevese 
etal. 1994, Cristiani et al. 1996, Di Clemente et al. 1996). 
As well as this connection with low-luminosity AGN, this 
effect would have important implications for the interpre- 
tation of the variability-luminosity relation: Neglecting the 
redshift/wavelength dependence would lead to an overes- 
timation of the real rest-frame optical variability of high- 2 
objects, introducing an artificial flattening in the variability- 
luminosity relation, given the strong coupling between L and 
z. 

In this paper we analyse the variability-luminosity- 
redshift relationship for the QSOs in the SGP sample. Our 
approach differs from that of Hook etal. (1994) in that we 
study variability as a function of both L and z using para- 
metric fits to the data. The advantage of model-fitting over 
a correlation analysis is that the uncertainties in the mea- 
sured quantities can easily be incorporated. Indeed, as we 
will demonstrate, photometric, sampling and k-correction 
uncertainties are a limiting factor in variability studies. 

Our main goals in this study are: 

• (i) to determine whether the data are consistent with 
a wavelength-dependent variability similar to that observed 
in low luminosity, nearby AGN, and 

• (ii) to test Poissonian (or 'Christmas-tree') models for 
QSO variability. In its simplest version this model predicts 
a 'l/v^/V' relationship between the rms variability and the 
mean luminosity. The proportionality constant carries infor- 
mation on the typical energy and time-scale of the individual 
pulses. More general models, in which the pulse properties 
scale with the QSO luminosity, are also investigated. 

Section 2 introduces the variability indices used in this 
work. Section 3 presents a quantitative discussion of the 
variability- wavelength relation in nearby AGN, setting a 
reference point for the results of the fitting analysis. The 
data set employed is discussed in Section 4, while Section 5 
presents the fitting method and a thorough discussion of the 
sources of error. The results of the analysis are presented in 
Section 6. The implications of our findings for Poissonian 
models are discussed in Section 7. Finally, in Section 8, we 
summarize our conclusions. 



2 VARIABILITY INDICES 

Two variability indices are used in this work. The first is 

„ = ^1 (1) 

Lb 

where Lb is the mean B-band luminosity in the light curve 
and <j(Lb) is its standard deviation. 

The rms variability is not an ideal index for samples 
spanning a wide range of redshifts. Because of the time dila- 
tion factor Atrest-frame = At bs/{1 + z), high-2 QSOs are 
observed for a substantially shorter rest-frame time than 
their low- 2 counterparts. Since QSO variability is correlated 
within time scales smaller than ~ 2-3 yr (Giallongo et al. 
1991, Hook etal. 1994, Cristiani etal. 1996), time dilation 



alone could lead to a systematic decrease of v for increasing 
z. 

One way to avoid this bias is to use an index based on 
observations separated by a fixed rest-frame interval, or by 
an interval longer than a critical value (see also Giallongo 
et al. 1991, Cristiani et al. 1996). With this purpose we define 
our second variability index as 



u = [L B (tj) -L_ B (U)f (2) 

Lb 

where t t and tj are epochs separated by more than 3 yr. 
This index corresponds to the asymptotic limit of the Struc- 
ture Function (SF) of each individual QSO. The SFs of indi- 
vidual SGP objects are of course poorly defined given that 
only seven epochs are available. The ensemble SF, however, 
rises steadily for intervals shorter than fa 2-3 yr, becoming 
essentially flat (though fairly noisy) for longer time-scales 
(Hook etal. 1994, Cristiani etal. 1996, Giallongo etal. 1991, 
Trevese et al. 1994). 



3 VARIABILITY VERSUS WAVELENGTH IN 
NEARBY AGN 

One of the goals of this paper is to compare the variability- 
wavelength relation of QSOs to that of low-luminosity, 
nearby AGN. In order to quantify this relationship we have 
parametrized it as a power law v(X) oc A - and used pub- 
lished data to estimate b for nearby objects. 

Among the five type-1 Seyferts studied by Edelson et al. 
(1990), four (NGC 4151, 3516, 5548 and Mrk 335) vary more 
at 1450 A than at 2885 A. We estimated b for these AGN 
from 

^ _ log (^145o/^2885) 

~ log (1450/2885) 

using «i45o and «2885 as listed in Table 1 of Edelson et al. 
(1990). The obtained values of b range from 0.47 to 0.66, 
i.e., the variations are between 40 and 60 per cent larger at 
1450 than at 2885 A. Interestingly, 3C273 (the only luminous 
AGN in their sample) also shows a wavelength-dependent 
variability with b fa 1. 

Aretxaga, Cid Fernandes & Terlevich (1996) estimated 
b from a comparison of the UV and optical light curves of 
NGC 4151 and 5548. They conclude that these nuclei vary 
2-3 times more in the UV (1300-1800 A) than in the B- 
band, which corresponds to b between 0.6 and 1. 

Finally, the results of Kinney et al. (1991) and Paltani & 
Courvoisier (1994) were also used to estimate b. This yielded 
0.6 < b < 2.9 (median b = 1.3) for Kinney etal.'s data 
(excluding blazars) and 0.2 < b < 1.7 for Paltani & Cour- 
voisier's data, but these values are less accurate than those 
quoted above. 

The observations therefore point to values of b between 
fa 0.5 and 1.5 for nearby AGN. While a v oc A - law is obvi- 
ously a simple tentative parametrization, the derived values 
of b provide an important reference point for the analysis of 
the v(X) relation in QSOs. 
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4 DATA: THE SGP SAMPLE 

The SGP sample consists of 283 QSOs monitored since 1975 
with the UK Schmidt telescope at the AAO (see Hook et al. 
1994 for details of the observations). Redshifts are in the 
0.43-4.01 range and apparent B magnitudes are between 
17 and 21. Eleven plates were obtained, defining 7 differ- 
ent epochs spanning 16 yr. The rest-frame light curves are 
between 3.2 and 14.5 yr long allowing for time dilation. Al- 
hough the sampling of each object is poor and the photomet- 
ric errors are not ideal (typically ±0.07 mag), the large num- 
ber of objects and wide coverage of the luminosity-redshift 
plane makes this sample one of the best QSO variability 
databases available. 

The data consist of B-band light curves and redshifts 
for all 283 QSOs. Six objects were excluded from the origi- 
nal list, five of which had unreliable (or missing) ras values. 
The other object is 0055-2659, a highly variable luminous 
QSO, which seems to be a peculiar source (see discussion in 
Hook etal., who also excluded 0055-2659 from their analy- 
sis). For each of the remaining 277 QSOs we computed Lb, 
v and u. The k-correction was computed as in Hook etal., 
except for the fact that we shall investigate spectral indices 
other than a = -0.5 (F v oc v a ). H = 50 kms _1 Mpc _1 and 
go = 0.5 are used throughout the paper, unless stated oth- 
erwise. Note, however, that v and u are independent of the 
k-correction and cosmological parameters, since they come 
in as multiplicative factors both in the numerator and de- 
nominator. Cosmology and k-correction affect only our es- 
timate of Lb, while v and u still measure the variability at 
A = Xb/(1 + z). 



5 FITTING METHOD AND ERROR 
ANALYSIS 

To examine the relationship between v, Lb and z we have 
used parametric fits in the (v, Lb, z) space. This is done by 
fitting the data with a function of the form 

v(L^,z) = ALB~ a (l + z) b (3) 

and analogously for u(Lb, z). The (l-\-z) b term is equivalent 
to (A/As) -6 , with A as the rest-frame wavelength in which 
the QSO is observed. The main advantage of equation (3), 
apart from its simplicity, is that its parameters can be used 
not only to test models for the v(L) relation (Section 7), but 
also to compare the v(X) relation of QSOs to that observed 
in lower luminosity AGN (Section 3). 

The fitting problem is linearized defining new variables: 

{/ = log v 
x = log (i s /10 12 i°) (4) 
y = log (i + z) 

Taking the logarithm of equation (3) we obtain 

f(x, y) = ax + by + c (5) 

We therefore seek to fit a plane to the data points in the 
space (x, y, /). The fit is performed using a \ 2 minimization 
algorithm which takes into account the errors in both / and 
x, but neglects errors in y. The \ 2 reads as 



2 _ (f< ~ aXt ~ ~ c ) 

i = l 

Any meaningful model-fitting requires realistic esti- 
mates of the errors involved. Several sources of error affect 
each data point in the sample. 

• Lb - the k-correction necessary to bring all the observed 
luminosities to a common reference wavelength (the B-band 
in this case) introduces a substantial uncertainty in Lb - The 
reason is that while a typical QSO spectrum needs to be as- 
sumed to compute the k-correction, individual QSOs have a 
wide range of spectral properties (e.g. Francis 1993). Red- 
dening may also affect Lb, although the indications are that 
QSOs suffer little intrinsic reddening (O'Brien, Gondhalekar 
& Wilson 1988, and references therein), and galactic redden- 
ing is likely to be low for the SGP objects. 

• v and u: the variability indices are affected by both 
the photometric errors and the sparse sampling of the light 
curves. 

• z: the broad width of the emission lines in QSOs, as 
well as possible peculiar motions, introduce uncertainties in 
the measured redshifts. These, however, are much smaller 
than the ones in v and Lb, and we shall neglect them. 

The uncertainties involved in deriving the mean lumi- 
nosity of QSOs are not often allowed for when analysing 
the variability-luminosity relation. Recognizing the presence 
of such uncertainties has a strong impact on the resulting 
slopes. This is discussed by La Franca et al. (1995), who show 
that underestimating or neglecting the errors in an indepen- 
dent variable (x in our case) leads to an underestimation of 
the true slope of the f(x) relation, i.e. to a small absolute 
value of a. 

In what follows we explain how we have estimated the 
uncertainties of Lb and v. 

5.1 Errors in the luminosity 

In the case of Lb we adopted the strategy of attributing the 
k-correction uncertainties to a range of possible spectral in- 
dices. Observed spectral indices in QSOs cover a wide range, 
from a PS —2 to +2, depending on which spectral window 
is used. The mean indices for the QSOs in the Large Bright 
QSO Survey (Morris etal. 1991, and references therein), for 
instance, are ~a = —0.5 for the 1455-2190 A window and 
— 0.8 for the 2190-3060 A window, with an rms dispersion 
<j a = 0.5 for both wavelength regions (Francis 1993). In this 
study we shall assume that a follows a Gaussian distribu- 
tion, centred at either a = —0.5 or —1.0, with a a = 0, 0.25 
or 0.5. For simplicity we follow the usual assumption that 
a is independent of both luminosity and redshift (but see 
Francis 1993). 

As the k-correction increases with redshift, the uncer- 
tainty in a propagates to an uncertainty in Lb (and x) which 
increases with z. The k-correction used also allows for Lya 
in emission as well as for absorption by intervening gas at 
A < 912 A (see Hook et al. 1994 for details). The uncertainty 
in x was thus computed with a small simulation, evaluating 
x for each QSO many times with a drawn from its Gaussian 
distribution. For a a = 0.25 the typical error-bars in x range 
from 0.08 to 0.17 for 1 < z < 4, corresponding to errors of 
18-39 per cent in Lb- 
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Figure 1. Extreme configurations for the errors in x (a) and 
log vq (b) . The luminosities in panel (a) were computed with ~ = 
-0.5. 

Each point in the light curve has a photometric error e m 
between 0.05 and 0.11 mag, depending on the object's appar- 
ent magnitude (see Table 2 in Hook et al. 1994). The mean 
luminosity obtained from the seven epochs is determined to 
an error \/7 better than the individual points. This com- 
ponent, although it was small compared to the k-correction 
uncertainty, was included the evaluation of e x perturbing 
each point in the light curve with a Ams drawn randomly 
from the true (apparent-magnitude dependent) error distri- 
bution of the SGP plates. The value of x was then computed 
1000 times for different realizations of the perturbations and 
spectral indices, yielding a reliable estimate of the combined 
photometric + k-correction uncertainties upon x. 

5.2 Errors in the variability index 

The photometric noise contribution to the rms variabil- 
ity was subtracted in quadrature from the measured rms: 
v = (v 2 - el f 12 , with e v = 2. 5 log e x e m . Six QSOs (all of 
which are in the high luminosity half of the sample) were ex- 
cluded from the analysis because they varied less than their 
estimated photometric error. 

The quadrature correction only improves our estimate 
of v; it does not provide a measure of how accurate this esti- 
mate is. The uncertainty in vo is most easily evaluated with 
a Monte Carlo simulation. For each object we: (1) perturb 
each point in the light curve by an error drawn from the ap- 
propriate error distribution, (2) compute v and vo, and (3) 
repeat the above steps 1000 times. The standard deviation 
of the resulting vo's gives a realistic estimate of e v o- 

The error estimates obtained above are entirely due to 
photometric noise. The poor sampling of the light curves 
ought to introduce a further uncertainty in vq. An estimate 
of the total photometric + sampling uncertainties upon vo 
can be obtained assuming that QSO variability is a Gaussian 



process, i.e. that is(t) is a Gaussian-random variable with 
mean Lb and standard deviation ao = vo x Lb- The proce- 
dure here was, for each object, to: (1) compute Lb, v and 
vo, (2) simulate a 7 points Gaussian light curve is(t) with 
standard deviation ao, (3) add a perturbation drawn from 
the appropriate error distribution and recompute vo and (4) 
repeat the above steps 1000 times, obtaining a new estimate 
for e„o. This new estimate measures the degree to which 
we can measure the correct standard deviation of a Gaus- 
sian process given only 7 points, all of which are affected 
by photometric errors. This method provides an upper limit 
to e„o, since QSO variability is not Gaussian in the time 
scales sampled by the SGP data, otherwise the ensemble SF 
should be flat, showing no correlated variability at all. On 
the contrary, the SF rises during the first 2-3 yr and only 
then flattens to a roughly constant level (Hook et al. 1994). 
Since many of the observations were carried out within (rest- 
frame) intervals smaller than 3 yr, we expect the observed 
fluctuations to be smaller than the ones corresponding to a 
truly Gaussian process. 

The values of e v o obtained with the 'Gaussian-method' 
are typically ~ 40 per cent larger than the ones obtained 
considering the photometric errors alone. The true uncer- 
tainties in e v o lie between these two estimates. The errors 
in / = log vo are easily evaluated either by using a stan- 
dard propagation of errors or by directly computing / in 
the simulations. 



6 RESULTS 

It should be clear from the discussion in Section 5.1 and 
Section 5.2 that though the sources of errors in luminosity 
(x) and variability index (/) are well understood, the exact 
magnitude of these errors is not well defined. It is there- 
fore important to investigate how sensitive the best fitting 
parameters a, b and c are to the choice of errors. We shall 
discuss nine different combinations of the errors in x and /: 

• e x : a a = 0, 0.25, 0.5 

• ef. Small, Large, Intermediate 

where 'Small', 'Large' and 'Intermediate' correspond to er- 
rors e v o computed using only the photometric errors, pho- 
tometric + sampling uncertainties, and the mid-point be- 
tween these two extremes, respectively. Fig. 1 illustrates 
some of these configurations. Less variable objects (low /) 
have systematically larger e/'s due to the increasingly uncer- 
tain quadrature correction to the variability indices. More 
luminous objects have larger e x , since they have higher red- 
shifts and are thus more affected by the k-correction uncer- 
tainties. In addition, two values for the mean spectral index 
a will be investigated: —0.5, the value used by Hook et al. 
(1994), and —1, the value used in Trevese et al. (1994). 

The data is shown in Fig. 2, where we plot / versus y, f 
versus x and y versus x. The left and right columns show the 
measurements with the smallest and largest possible errors 
allowed in this study, respectively. Fig. 3 attempts a visu- 
alization of the variability in the luminosity-redshift plane. 
Each QSO is represented by a square whose size scales with 
vq. The anticorrelation between variability and luminosity 
is clearly visible: squares are smaller towards larger x. The 
variability-redshift relation shows a large scatter, though the 
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Figure 2. log-uo X y, log-uo X x and x X y for the QSOs in the 
SGP sample. Error bars correspond to (•,.••) = (0, Small) for 
the plots in the left and (0.5, Large) for those in the right. 

tendency is for vo to increase with y within fixed luminosity 
bands. 

The results of the parametric fits are presented in Ta- 
ble 1. The best fitting parameters a, b and c are very sensitive 
to the choice of errors. In general a gets increasingly nega- 
tive and b increases as a a (and consequently e x ) increases, 
whereas the opposite happens as t/ increases. The param- 
eters a and b always vary in opposite directions due to the 
strong coupling between luminosity and redshift. The inde- 
pendent term c tends to follow a in its dependence on e x 
and Cf. While a and c are nearly unaffected by the mean 
spectral index used in the k-correction, b increases as a goes 
from —0.5 to —1. 

The changes in the best-fitting parameters resulting 
from different choices of errors are substantial. For a = —0.5, 
for instance, the pair (a, b) changes from ( — 0.81,1.52) to 
( — 0.19, 0.12) as we move from the case in which e x is large 
and Cf is small to the opposite extreme. The errors in the 
fit parameters (see Fig. 4) are smaller than this range. It is 
therefore clear that the uncertainty in a, b and c is domi- 
nated by our lack of detailed knowledge of the magnitude of 
the errors involved. On the assumption that Table 1 covers 
the whole spectrum of possibilities for the errors, and dis- 
carding cases where either / or x carries essentially all the 
errors (i.e., maximum e/ and minimum e x and vice- versa), 
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Figure 3. Projected representation of the variability in the x X y 
plane. The size of each square scales linearly with vq . The six 
objects marked with filled squares are the ones which vary less 
than the expected photometric error, with vq represented by \v 2 — 
e^l 1 / 2 . The absolute magnitude corresponding to Lg is shown in 
the upper axis. 

we can safely bracket the best-fitting parameters within the 
following intervals: 

-0.6 < a < -0.2 
0.1(0.2) < b < 1.0(1.3) (6) 
-1.0 <c< -0.7 

where values of b outside parentheses correspond to a = 
— 0.5 fits, and those inside are for a = —1. Unless we choose 
a specific error-configuration, it is not possible to determine 
a, b and c to a better degree of precision. 

Despite the large ranges in the fitting parameters, two 
strong conclusions can be derived from our analysis. First, 
we confirm that variability and luminosity are anticorrelated 
in the SGP QSOs (a < in all fits), as previously found by 
Hook et al. (1994 — see also Pica & Smith 1983, Cristiani, 
Vio & Andreani 1991, Cristiani et al. 1996). Secondly, the 
SGP data point to the existence of a positive correlation be- 
tween variability and redshift (b > in all fits), in agreement 
with Cristiani et al. (1996), who used a different method of 
analysis. Furthermore, the obtained values of b are consis- 
tent with the ones estimated from the variability-wavelength 
relation observed in less luminous AGN (0.5 < b < 1.5). 
This is a strong indication that the variability-wavelength 
anticorrelation observed in nearby AGN is also followed by 
high-redshift QSOs (see also Paltani & Courvoisier 1994, Di 
Clemente et al. 1996). 

The uncertainties in a, b and c are illustrated in Fig. 4, 
where we plot the one, two and three sigma confidence levels 
in the a versus b and a versus c planes for the four ex- 
treme combinations of a a and ef. (a a ,Cf) = (0, Small), 
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-1.0 


0.50 


S 


-0.78 


1.87 


-1.14 


1.751 


-1.0 


0.50 


I 


-0.59 


1.32 


-1.03 


1.464 


-1.0 


0.50 


L 


-0.46 


1.00 


-0.98 


1.176 



Table 1. Best fitting parameters a, b and c and reduced x 2 f° r 
the log^o = ax + by + c fits. The different lines correspond to 
different combinations of ~, <7 a and t^. S, 1 and L denote small, 
intermediate and large errors in / = log^o. 
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Figure 4. Ax 2 maps for the v(Lb, z) fitting parameters, showing 
the 68.3, 95.4 and 99.7 per cent confidence level contours. Num- 
bers in parentheses are {-.',. tf). All contours are for ~ = —0.5. 

(0.5, Small), (0, Large) and (0.5, Large). Other combina- 
tions of the errors in x and / yield contours within the gen- 
eral area delimited by these cases. The goodness of the fits, 
measured by the reduced % 2 , are listed in Table 1. An ac- 
ceptable x 2 here means x 2 /rf.o./ < 1.288 given the 268 de- 
grees of freedom — the probability of larger values occurring 
by chance is 10 -3 . Only two of the fits produce acceptable 



X 2 values, notably (but not surprisingly) those with large e x 
and tf. In any case, though the uncertainties in a, b and c 
are large, it is clear that a > and b < models are ruled 
out by the observations. 

6.1 Tests of the fits 

A few variations over the fitting analysis were tried to test 
the sensitivity of the fits to some aspects of the data. 

The k-correction includes a correction for absorption 
by the LyajLyfi forest and for Lyman-limit systems (Irwin, 
McMahon & Hazard 1991). Since this correction affects only 
the QSOs with z > 3.28, we have experimented with exclud- 
ing the 17 such objects from the sample to test how strongly 
they affect the fits. The resulting changes in a, b and c were 
negligible, typically less than 0.1. 

Fits were also made to subsets of the data dividing the 
sample into luminosity and redshift halves. The values of 
a for the lower luminosity and redshift half-samples were 
systematically less negative than the values for the high lu- 
minosity and redshift halves. As usual, b showed the oppo- 
site behaviour. Whilst this result indicates a bending of the 
f(x,y) plane, the major factor affecting the fits is the in- 
crease in the errors as we go towards more luminous, high 
redshift objects. In fact, the low-i and low-2 fits are closer 
to the ones obtained for the full sample, since these objects 
weight more in the fits as they have smaller errors. A wider 
coverage of the luminosity-redshift plane and smaller errors 
would be required to properly assess the reality of this bend- 
ing. 

We also computed fits neglecting the presence of errors 
in x and /. These yielded a = —0.26 and c = —0.87, whereas 
b was 0.19 for a = —0.5 and 0.31 for a = —1. Finally, 
fits were also performed for Ho = 100 kms -1 Mpc -1 . As 
expected, the coefficients a and b are not affected by Ho, 
but c is systematically smaller by Ac = a log 4 than for 
Ho = 50 kms -1 Mpc -1 . This difference amounts to PS 0.1- 
0.5 for the values of a in Table 1. 

6.2 Fits with a fixed wavelength/redshift 
dependence 

Instead of allowing for a wavelength dependence in the fits, 
one could alternatively start with the assumption that QSOs 
follow the same variability-wavelength relationship observed 
in nearby AGN. Adopting a fixed value of b, we can apply a 
correction to the observed variability indices to convert them 
to a common wavelength for all QSOs and then examine the 
consequences for the variability-luminosity relationship. 

In Fig. 5 we plot the best-fitting values of a and c as a 
function of b for different combinations of a a and tf. The 
tendency is for both a and c to decrease with increasing b. 
Adopting a value of b similar to those found in nearby AGN 
(0.5 < b < 1.5) would result in -0.8 < a < -0.3, while c 
would lie in the —1.2 to —0.8 interval. The A% 2 profiles are 
also shown in Fig. 5 to illustrate the change in the quality 
of the fit for different values of b. 

Note that even for 6 = 0, i.e., no wavelength depen- 
dence, the slope a is critically dependent on the choice of er- 
rors, particular for the luminosity. The change in a between 
fits with a a = and 0.5 is comparable to that resulting from 
changing b from to 1. 
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Figure 5. Best-fitting parameters a and c as a function of b for x = lo § ( l b/ 10 ' 3 l b) 

fixed 6 fits for different combinations of the errors in luminosity 

O) and variability index (/). The top panel shows the change in Figure 6. As Fig. 3, but for variability index u a . 

the x 2 with respect to the model with minimum x 2 • The range 
of values of b found in nearby AGN is indicated for comparison. 



6.3 Results for the SF-based variability index 

Fits analogous to those for v were also performed for the SF- 
based variability index, u, defined in equation (2). The anal- 
ysis of the errors in u was carried out as for v, with the excep- 
tion that the quadrature correction to u is Mo = (m 2 — f 2 ) 1 ^ 2 , 
with t u = 5 log e x e m . This time, 10 objects were excluded 
from the analysis for having u < t u . These are marked with 
filled squares in Fig. 6, where, analogously to Fig. 3, we rep- 
resent the variability in the x x y plane with squares whose 
sizes scale with Mo. We find that e„o is typically 10-20 per 
cent larger than e„o, and that, as for vq (Fig. 1), objects 
with small Mo are subjected to larger errors. Fig. 7 shows 
the variability versus luminosity and redshift relations for 
Mo in the cases of minimum and maximum errors. 

The results of the fits are listed in Table 2, and con- 
fidence contours are plotted in Fig. 8. It is noticeable that 
a tends to be more negative and b more positive than for 
the Va-Hts. This happens essentially because log Mo covers 
a larger dynamic range than does log^o, which naturally 
implies larger slopes. The range of log Mo in the sample is 
from —2.1 to —0.1, with an rms dispersion of 0.3, whereas 
for log Vq the range is —1.6 to —0.2 and the rms is 0.2. No- 
tice, however, that this difference is only noticeable for fits 
with large a a . Fits with no errors resulted in a = —0.24 and 
c = -0.80, with b = 0.05 for a = -0.5 and 0.17 for a = -1. 

Fig. 9 shows the ratio between the two variability in- 
dices as a function of no, z and x. The index Mo is typically 
~ 20 per cent larger than vo since it is based on variations 
over a longer time base-line, but there is no obvious trend of 
the ratio uq/vq with vo, z or x. This indicates that vo is not 
severely affected by the time-dilation bias discussed in Sec- 
tion 2, otherwise uq/vq should exhibit a systematic increase 
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Co 


e f 


a 


b 


c 


X 2 /d.o.f. 


-0.5 


0.00 


S 


-0.26 


0.27 


-0.71 


3.980 


-0.5 


0.00 


I 


-0.19 


0.08 


-0.73 


2.375 


-0.5 


0.00 


L 


-0.18 


0.04 


-0.75 


1.771 


-0.5 


0.25 


S 


-0.48 


0.77 


-0.86 


3.430 


-0.5 


0.25 


I 


-0.30 


0.33 


-0.79 


2.320 


-0.5 


0.25 


L 


-0.25 


0.19 


-0.78 


1.749 


-0.5 


0.50 


S 


-1.13 


2.23 


-1.30 


2.077 


-0.5 


0.50 


I 


-0.81 


1.48 


-1.11 


1.814 


-0.5 


0.50 


L 


-0.61 


1.02 


-1.01 


1.528 


-1.0 


0.00 


S 


-0.25 


0.36 


-0.70 


3.946 


-1.0 


0.00 


I 


-0.20 


0.17 


-0.72 


2.385 


-1.0 


0.00 


L 


-0.19 


0.17 


-0.75 


1.769 


-1.0 


0.25 


S 


-0.50 


1.04 


-0.86 


3.447 


-1.0 


0.25 


I 


-0.30 


0.48 


-0.79 


2.340 


-1.0 


0.25 


L 


-0.24 


0.32 


-0.79 


1.741 


-1.0 


0.50 


S 


-1.12 


2.76 


-1.29 


2.077 


-1.0 


0.50 


I 


-0.82 


1.92 


-1.12 


1.830 


-1.0 


0.50 


L 


-0.61 


1.33 


-1.01 


1.523 



Table 2. As Table 1 but for the logiio = ! - ' ' * * fits. 

with z. This is not seen in Fig. 9, with the possible exception 
of z > 3 objects. This, however, is not necessarily a general 
result, as it could be the result of the particular distribution 
of time intervals in the SGP sample. The scatter in Mo /vo 
increases substantially for luminous, high-2, less variable ob- 
jects, most probably due to the larger errors, particularly in 
Mo. Because of time-dilation, the number of pairs of plates 
involved in the computation of Mo decreases from 15 to 7 
as z increases from 0.4 to 3.6, and for the highest redshift 
object in the sample Mo is computed from only two epochs. 
Analogously to equation (6), we can bracket the best- 
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a= 0.0 e, = Small u= 0.5 e, = Large 




Figure 7. As Fig. 2 but for uq as the variability index. 




Figure 8. As Fig. 4, but for the uo(Lg,z) fits. 



fitting parameters for the Mo fits between the following in- 
tervals: 



-0.8 < a < -0.2 
0.1(0.2) < b < 1.5(1.9) 
-1.1 <c< -0.7 



(7) 



(where values of b outside parentheses correspond to a = 
— 0.5 fits and those inside are for a = —1.0). Note, however, 
that the \ 2 values in Table 2 are larger than those for the 
«o-fits, which indicates that our fitting equation is too sim- 
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Figure 9. Ratio of the two variability indices, uq and vq as a 
function of vq , redshift and luminosity. Filled circles indicate ob- 
jects for which v < e v and/or u < e u . The upper scale in the top 
panel shows the absolute magnitude corresponding to x. 



plistic and/or other quantities besides L and z are involved 
in defining the variability properties of QSOs. 



7 TESTING MODELS FOR QSO VARIABILITY 

So far we have focused on the analysis of the data. In what 
follows we shall demonstrate that, despite the uncertainties 
discussed, the SGP sample imposes interesting constraints 
on theoretical models for AGN variability. 

Testing the consistency of Poissonian models for the op- 
tical variability of QSOs is of major interest in this study. 
In this general class of models variability is regarded as 
beibg the result of the random superposition of uncorre- 
cted pulses. In this section we use the results of our fits to 
test the validity of this scenario and to constrain its basic 
parameters. 

The essential elements of a Poissonian model are the 
rate, energy and time-scale of the pulses. These are denoted 
by v, Eg and r B respectively, the subindex indicating that 
we are referring to the B-band. The hypothesis that the 
events are randomly distributed in time yields the following 
expressions for the mean B-band luminosity and net rms 
deviation: 



Lb = vE P B l ' f V ar 

and 

vb = flLl (Lb/L p b 



-1/2 



(8) 



(9) 



where L g = E g /r g denotes the typical luminosity of an 
individual pulse, with E g and r g as the mean B-band en- 
ergy and lifetime of the pulses respectively (Cid Fernandes 
1995). Equations (8) and (9) allow for the possibility that a 
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fraction fbck = 1 — fvar of the mean luminosity Lb is due 
to a constant 'background' component in the light curve. It 
is clear that equation (9) is simply a 'one over square root 
of A' law, diluted by the fact that a fraction fbck of the 
observed luminosity is not variable. In what follows we use 
the SGP data to constrain fbck, E P B and v. 
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7.1 The background fraction: fbck 

Several sources can act as a constant background in 
AGN light curves: starlight from the host galaxy, a nu- 
clear/circumnuclear star cluster, scattering, and/or the non- 
variable part of the nuclear engine. Starlight from the host 
galaxy is certainly present in the light curves of QSOs, but 
its contribution to the total luminosity is likely to be small, 
as indicated both by direct imaging studies (Aretxaga, Boyle 
& Terlevich 1995, Hutchings 1995) and by the absence of 
strong stellar features in the spectra of QSOs. The nuclear 
engine itself, however, could well have a non-variable com- 
ponent. This is, for instance, the case in the starburst model 
for AGN; in this model a substantial fraction of the optical- 
UV luminosity is the result of the presence of young stars 
(Aretxaga & Terlevich 1994). In an accretion disc model, 
one can also expect that part of the disc will remain un- 
affected by the instabilities promoting the variability, thus 
playing the role of a background component. Scattering of 
the nuclear photons in an extended region, a central ingre- 
dient in current unified models (e.g. Antonucci 1993), would 
smooth out the variations and would also appear as a con- 
stant background. 

An upper limit for the background fraction can be ob- 
tained directly from the data using the average and min- 
imum luminosities in the light curves. Since Lg' n > L b g k 
and L b g k = fbckLg, we have that 



fbck < — — 



Fig. 10 shows the results of this exercise. As expected, 
Lg' n I Lg increases with luminosity, since the chances of a 
light curve reaching its true minimum decrease as the rate 
of events increases. Accordingly, the less luminous objects 
are those with a smaller upper limit for fbck- Taking this 
systematic effect into account, we estimate that fbck < 0.7 
for all QSOs. Note that several points in Fig. 10 lie above 
this limit. This conservative estimate allows for the possi- 
bility that Lg' n , which is evaluated from a single point in 
the light curve, may be affected by the photometric errors. 
This is an important model-independent quantity for vari- 
ability studies. If QSOs do contain a constant light source 
aside from the variable component, then this non- variable 
source contributes at most 70 per cent of their B-band lu- 
minosity. As a comparison, the ratio of minimum to average 
B luminosities in NGC 4151 and 5548 is ftj 0.5 in both cases, 
implying fbck < 0.5 (Aretxaga & Terlevich 1993, 1994, Cid 
Fernandes, Terlevich & Aretxaga 1996). Given the lower rate 
of events in these objects, this estimate is probably closer to 
its true value for them than it is for QSOs. 
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Figure 10. Upper limits for the background fraction (fbck) f° r 
the QSOs in the SGP sample. 



7.2 The Simple Poissonian model: Constraints on 
the energy of the pulses 

The simplest version of a Poissonian model is one in which 
the energies and time-scales of the pulses and the back- 
ground fraction are the same for all objects, regardless of 
their luminosities, redshifts, etc. We shall refer to this case as 
the Simple Poissonian model. This model makes two strong 
testable predictions: (1) the slope of the variability luminos- 
ity relation is a = —1/2, and (2) the energies and time-scales 
of the pulses are similar for objects from Seyfert galaxies to 
QSOs. 

Our analysis of the SGP data showed that a is depen- 
dent on the rather uncertain magnitude of the errors in L 
and, to a lesser extent, on the variability index. Interest- 
ingly, the value —0.5 is nearly at the centre of the range of 
values resulting from different combinations of errors. We 
therefore regard these observations as being consistent with 
the Simple Poissonian model, though it is clear that we are 
not able to test definitively the a = —0.5 prediction with 
the data presently available. 

On the assumption that the model is valid, the inde- 
pendent term c can be used to estimate the typical energy 
of the pulses. Taking the logarithm of equation (9) 



log VI 



^log f var 



and comparing it to / = ax + by + c computed for y = 
(i.e. z = and A = As), we derive 



1 / L p 

'" = ^ l°g fva B 



2 



10 12 i® 



(10) 



The v(Lb,z) fits of the SGP sample yielded values of 
c between —1.1 and —0.8 for those fits where a ftj —0.5, 
whilst the fits for the SF-based index u(Lb, z) resulted in a 
value of c between —1.1 and —0.7, also for fits yielding a ftj 
— 0.5. The theoretical prediction (equation 9) applies to the 
rms variability as measured from independent, uncorrelated 
points in the light curve, so in what follows we use the latter 
constraints for c. We then have that 



6.3 x 10 9 ig < f var L p B < 4 x 10 10 ig 



(11) 



These limits, coupled with the constraint that 0.3 < 
fvar < 1 yield limits on the pulse luminosity of 



6.4 x 10 9 ig < L P B < 1.3 x 10 n ig 



(12) 



10 R. Cid Fernandes, I. Aretxaga, R. Terlevich 



The value of c does not by itself provide independent 
constraints on the average energy and time scale of the 
pulses. We can nevertheless use the time-scales derived from 
the SF analysis in Hook et al. (1994). Their study has shown 
that the ensemble SF of the QSOs in the SGP sample rises 
for At < 2-3 yr and flattens off for larger lags, indicating that 
QSO variability is uncorrelated on longer time-scales (see 
also Cristiani et al. 1996). The pulse lifetime must therefore 
be in the neighbourhood of 2 yr, since the SF of a light curve 
resulting from the random superposition of many pulses is 
equal to the mean SF of the individual pulses. Assuming 
that the pulse lifetime is in the 1.5 < r g < 3 yr range, then 
yields the following constraints on the pulse energy: 

1.5 x 10 50 erg <~E% < 6.3 x 10 51 erg (13) 

This range incorporates the constraints in c, f var and 
r B , all of which are based on conservative estimates. The re- 
sulting limits on the pulse energy can therefore be regarded 
as a robust constraint. These limits are nevertheless cos- 
mology dependent: Ho = 100 kms -1 Mpc -1 would result in 
energy limits 4 times smaller than those in (13). 
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a 


c 


n 


v /yr 1 


erg 


0.5 


-0.25 


-0.75 


61 


7.6 


3.8 


0.0 


-0.50 


-0.90 


231 


28.9 


1.0 


-0.5 


-0.75 


-1.10 


1112 


139.0 


0.2 



Table 3. Estimated parameters for generalized Poissonian models 
with 7 = 0.5, and —0.5. Pulse rates and energies correspond to 
a Mb — —26 QSO, and assume a background fraction of 50 per 
cent and a pulse lifetime of 2 yr. 

estimate typical pulse energies for the three models consid- 
ered, let us take f var = 0.5 and r g = 2 yr as fiducial values 
for an Lo QSO. The resulting pulse energies are listed in 
Table 3, along with the required pulse rates, derived from 
v = Lo fvar/ Eg. Note that the pulse energy for a Simple 
Poissonian model (7 = 0) is consistent with (13), as ex- 
pected. The uncertainty in c is ±0.1, which propagates to 
an uncertainty of 46 per cent for no, v and E g . 



7.3 Generalized Poissonian models 

The simple model, while attractive and consistent with the 
data presently available, is not the only possible Poissonian 
model for QSO variability. More generally, one can envisage 
a situation in which the pulse properties (energy and time- 
scale) depend upon the global luminosity of the object. In 

this case the vb oc Lb ^ prediction would no longer be 
valid. 

Consider a case where A = f V arL g depends on the lu- 
minosity of the object as A(Lb) oc Lb' 1 , the exponent 7 
being determined by the physical process which links the 
properties of the individual pulses to those of the object as 
a whole. The variability luminosity relation then becomes 

(7 — 1 ) /2 

vb oc Lb ■ The fit coefficient c can now be used to 

calibrate A(Lb), analogously to what was done for the Sim- 
ple Poissonian model (7 = 0). Let us use as a reference 
point the luminosity Lo of a 'typical' QSO, corresponding 
to Mb = —26, and write 

A(I ^ ) = Ao (if) 7 (14) 

We shall restrict our discussion to three possible values 
of 7: 0.5, and —0.5. These models would result in v(L) 
slopes a = —0.25, —0.5 and —0.75, respectively. While the 
data analysis showed that we cannot as yet decide which 
of these possibilities is to be favoured, we can safely reject 
values of a outside this range. For each value of a we now 
estimate the corresponding value of the independent term 
c from Fig. 8 and compute Ao. The results are listed in 
Table 3, where instead of Ao we list no = io/Ao, i.e., the 
number of individual pulses which have to be added up to 
account for the luminosity of a Mb = —26 QSO. The actual 
mean number of pulses present at any moment is Ao = 
fvar n o, which reduces to no in the absence of a background 
component. 

The values of no obtained give us information on the 
product fvarE g /r g at the calibration point Lo. In order to 



7.4 Discussion 

The most appealing aspect of the results presented above is 
their generality. The only assumption involved is that the 
pulses occur randomly in the light curve. While this leaves 
out cyclical/periodic processes, the analysis encompasses a 
large class of possible physical models. 

Gravitational microlensing has recently been explored 
as an external source of variability in QSOs (Hawkins 1996). 
As the lensing events are probably independent, the whole 
process can be seen as Poissonian, though the luminosity of 
the QSO itself would act as a strong background. While this 
scenario predicts achromatic variability if all the emission 
arises from regions of the same scale-size, some chromati- 
cism is expected if the continuum at differing wavelengths 
originates from different parts of an accretion disc (Wams- 
bsganss & Paczyhski 1991, Alexander 1995). However, the 
size of the Broad Line Region in the canonical picture is too 
large to be magnified by the microlensing action of solar- and 
sub-solar mass stars at cosmological distances (Schneider & 
Wambsganss 1990), and no flux variations are expected in 
the emission lines. While this model is not ruled out as a 
general mechanism of variability in high-redshift QSOs (but 
see Alexander 1995), the variability properties of the major- 
ity of low-redshift QSOs clearly depart from the predictions 
of this scenario. Of the six well studied variable (z < 0.3) 
PG QSOs, four show line flux variations (Maoz et al. 1994). 
More generally, the higher variability of low-redshift, low- 
luminosity QSOs is contrary to the microlensing hypothe- 
sis, as the optical depth of microlensing objects would not 
be sufficient to produce the observed fluctuations (Wambs- 
ganss 1990, Lewis & Williams 1996). 

One well-studied example of a Simple Poissonian model 
applied to AGN is the starburst model (Terlevich et al. 
1992, Aretxaga & Terlevich 1994, Aretxaga et al. 1996). This 
model makes predictions for energies (3-6 x 10 50 erg in the 
B band) and time-scales (1-4 yr) of the pulses, as well as 
for the background luminosity (fbck ~ 40 per cent). The su- 
pernova rate goes from PS 3 yr -1 for an Mb = —23 QSO to 
200 yr -1 for a QSO with Mb = —28. These predictions are 
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thoroughly consistent with the constraints derived in Sec- 
tion 7.2 (we refer the reader to Aretxaga & Terlevich 1994 
and Aretxaga et al. 1996 for a discussion of this model). In 
the context of the super-massive black-hole scenario, AGN 
variability is traditionally associated with accretion instabil- 
ities (e.g. Rees 1984, Abramowicz 1991, Wallinder, Kato & 
Abramowicz 1992). If random accretion events are respon- 
sible for the variations observed in QSOs then they must 
produce B-band energies of the order of 10 51 erg. For an 
efficiency of 10 per cent this yields a mass of 0.05 Mq if 10 
per cent of the photons produced come out in the B-band. 
Furthermore, about 30 such events must happen per year 
in an Mb = —26 QSO, so the accretion rate in the form of 
inhomogeneities must be about 1.5 Mq yr _1 . 

Whilst Table 3 provides useful order of magnitude esti- 
mates of the rate and energy of the pulses in an Mb = —26 
QSO, we cannot extrapolate them to other luminosities 
without detailed knowledge of individual dependences of 
fvar, Eg and r g on the total luminosity. Such information 
would be very important for distinguishing between differ- 
ent scenarios. If, for instance, variability is associated with 
accretion onto a black-hole, one might perhaps associate the 
pulse duration with an orbital or viscous time-scale. If the 
accreted objects always have the same mass, one might ex- 
pect the energy of the pulses to scale with the black-hole 
mass. Both time-scale and energy would then change as a 
function of luminosity. 

Determining how the pulse properties scale with the 
object properties should therefore be a major goal of future 
variability studies. In principle, the pulse lifetime can be em- 
pirically determined from the SF of individual objects, but 
a precise measurement requires long and well sampled light 
curves. The evidence so far is that the time-scale as derived 
from the SF does not vary strongly as a function of lumi- 
nosity (Hook et al. 1994, Cristiani et al. 1996). Pulse energies 
can only be directly measured if the events do not overlap. 
While this exercise can be tried in low luminosity AGN, in 
QSOs a large number of pulses are present at any time. This 
can be seen in Table 3: computing No ~ vt^ yields at least 
15 coexisting events in an Mb = —26 QSO. In such cases 
the pulse energy has to be derived indirectly, as above. Sim- 
ilarly, the background fraction cannot be determined when 
overlapping occurs, though it may be constrained through 
the method discussed in Section 7.1. 



8 SUMMARY 

The South Galactic Pole database was used to investigate 
the relationship between variability, luminosity and redshift 
in QSOs. This was done using parametric fits to the data in 
the 3D space defined by these three variables. Two vari- 
ability indices were used: v, the net rms variability, and 
u, an index based on observations carried out more than 
3 yr apart. This latter index avoids the problems related to 
time-dilation effects, but is subjected to larger errors than 
v. The variability-luminosity-redshift relation was fit with a 
v oc L a (l + z) function, and similarly for u. This approach 
to the problem allows a careful examination of the effects 
of the errors involved. The main results of the data analysis 
can be summarized as follows: 



• (i) We confirm Hook et al. 's (1994) result that vari- 
ability is anticorrelated with luminosity in the SGP QSOs 
(a < in all fits). 

• (ii) We also found that variability increases with red- 
shift (b > 0). This is in agreement with the results of Cris- 
tiani et al. (1996), but not with Hook et al.'s conclusion that 
the redshift correlation is a by-product of the coupling of 
L and z. In particular, we found that the values of b ob- 
tained for the SGP sample are consistent with those derived 
from the variability-wavelength anticorrelation observed in 
nearby AGN, indicating that this anticorrelation extends to 
luminous, high-2 objects (as first suggested by Giallongo 
etal. 1991). 

• (iii) The logarithmic slope of the v x L relationship was 
found to lie in the a = —0.7 to —0.2 interval, depending 
on the magnitude of the errors involved. Similar slopes were 
obtained for the SF-based index u(L). The uncertainties as- 
sociated with sampling and particularly the k-correction do 
not allow a more precise determination of this slope. 

The results of the data analysis were compared to theo- 
retical predictions for Poissonian models, understood as sys- 
tems in which the light curve is the result of a random su- 
perposition of pulses plus a possible background component. 
This comparison produced the following constraints which 
must be satisfied by physical models for the origin of the 
pulses: 

• (i) The minima in the light curves were used to estimate 
the contribution of a constant 'background' source to the 
total B-band luminosity of the QSOs, yielding fbck < 70 per 
cent. 

• (ii) We explored the hypothesis that QSO variability 
can be explained in a Simple Poissonian model, in which 
the energy and time-scales of the pulses are the same for all 
QSOs independent of their luminosity. This model predicts 
a 'one over square root of TV' relation between the net rms 
variability and the luminosity. We were able to put hard 
limits on the typical B-band luminosity of the pulses. These 
can then be coupled with the 1.5 < Tg < 3 yr limits on the 
time-scale of the pulses, as derived from the analysis of the 
ensemble SF, to constrain the B-band energy of the pulses 
between 1.5 x 10 50 and 6.3 x 10 51 erg. 

• (iii) Generalized Poissonian models, in which the pulse 
properties scale with the QSO luminosity, were also investi- 
gated. Models in which the pulse and object luminosities are 
related by L V B oc Lb 1 are also consistent with the SGP data 
provided -0.5 < 7 < 0.5. For an Mb = -26 QSO the pulse 
rates and B-band energies required in such models were all 
found to be within an order of magnitude of 30 yr -1 and 
10 51 erg respectively. 

Future observations should concentrate on reducing the 
errors both in variability and in luminosity, and on improv- 
ing the coverage of the L-z plane. Spectrophotometry and/or 
multi-colour monitoring of QSOs could in principle provide 
more stringent tests of the wavelength dependence of vari- 
ability. This would automatically remove the wavelength ef- 
fect, something which our simple v(X) oc X~ b law could only 
do in an approximate way. Residual correlations with red- 
shift, if any, would then have to be explained in another way 
(for instance, as due to cosmic age). It is also important to 
encourage the continuation of monitoring programs such as 
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that of Hook et al. (1994), since their wide time baseline 
provides the type of long term information which we cannot 
hope to obtain in the near future. Also, as the length of the 
light curves increases the errors associated with sampling 
effects should decrease. It would also be interesting to study 
the dependence of variability on properties other than L and 
z. For instance, radio-loud objects seem to vary more than 
their radio-quiet counterparts (Pica & Smith 1983, Maoz 
et al. 1994). Radio observations of the SGP objects would 
be important to verify this property and to better study the 
v-L-z relationship in these two families of AGN. 

Until now variability studies of high- and low-luminosity 
AGN have been largely treated as separate topics. In order 
to further our understanding of AGN, it will be of great im- 
portance to link the variability properties of QSOs to those 
of Seyferts and LINERs. Only then will we be able to address 
fundamental questions such as: Do low luminosity AGN fol- 
low the same variability-luminosity-wavelength relation as 
QSOs? Is the process Poissonian? What is the rate, energy 
and time-scales of the pulses and how do they vary from 
object to object? Answering these questions would narrow 
considerably the range of physical models for AGN variabil- 
ity. 

Acknowledgments: We are greatly indebted to Isobel 
Hook, Richard McMahon, Mike Irwin and Brian Boyle for 
making the SGP data available to us. Discussions with Mike 
Irwin, Geraint Lewis, Richard McMahon, Onno Pols and 
Laerte Sodre are also acknowledged. The hospitality of ESO 
(Santiago), where this work was completed, was much appre- 
ciated. The work of RCF was supported by CAPES/Brazil 
under grant 417/90-5. IA's work is supported by the EEC 
HCM fellowship ERBCHBICT941023. 



Hook I. M., McMahon R. G., Boyle B. J., Irwin M. J 1994, MN- 
RAS, 268, 305 

Irwin M., McMahon R. G., Hazard C. 1991, in The Space Distri- 
bution of Quasars, Crampton D. ed, ASP Conf. Ser. Vol. 21, 
(San Francisco), p. 117 

Kinney A. L., Bohlin R. C, Blades J. C, York D. G 1991, ApJS, 
75, 645 

La Franca F., Franceschini A., Cristiani S., Vio R. 1995, A&A 
299, 19 

Lewis G. F. & Williams L., 1996, preprint 

Maoz D., Smith P. S., Jannuzzi B. T., Kaspi S., Netzer H. 1994, 
ApJ, 421, 34 

Morris S. L., WeymannR. J., Anderson S. F., Hewett P. C, Foltz 

C, Chaffe F. H., Francis P. J., MacAlpine G. M. 1991, AJ, 
102, 1627 

Netzer H., Sheffer Y. 1983, MNRAS, 203, 935 

O'Brien P., Gondhalekar P., Wilson R. 1988, MNRAS, 233, 801 

Paltani S., Courvoisier T. J.-L. 1994, A&A, 291, 74 

Pica A. J., Smith A. G 1983, ApJ, 272, 11 

Rees M. 1984, ARAA, 22, 471 

Schneider P. & Wambsganss, J., 1990, ApJ, 237,42. 
Smith P. S., Balonek T. J., Elston R., Heckert P. A. 1987, ApJS, 
64,459 

Terlevich R., Tenorio-Tagle G., Franco J., Melnick J. 1992, MN- 
RAS, 255, 713 

Trevese D., Pitella G., Kron R. G., Bershady M. A. 1989, AJ, 98, 
108 

Trevese D., Kron R. G., Majewski S. R., Bershady M. A., Koo 

D. C. 1994, ApJ, 433, 494 

Uomoto A. K., Wills B. J., Wills D. 1976, AJ, 81, 905 
Wallinder F., Kato S. & Abramowicz M. A. 1992, A&AR, 4, 79 
Wambsganss, J., 1990, Ph.D. Thesis, Miinchen. 
Wambsganss, J. & Paczyhski B. 1991, AJ, 102, 864. 



REFERENCES 

Abramowicz M. A. 1991, in Duschl W. J., Wagner S. J., Camezind 
M., eds, Lecture Notes in Physics, 377, Variability of Active 
Galaxies, Springer- Verlag, Berlin p. 255 
Alexander T. 1995, MNRAS 274, 909 
Angione R. J. 1973, AJ, 78, 353 
Antonucci R. 1993, ARA&A, 31, 473 
Aretxaga I., Terlevich R. 1993, A&SpSci, 205, 69 
Aretxaga I., Terlevich R. 1994, MNRAS, 269, 462 
Aretxaga I., Boyle B., Terlevich R. 1995, MNRAS, 275, L27 
Aretxaga I., Cid Fernandes R., Terlevich R. 1996, MNRAS sub- 
mitted 

Bonoli F., Braccesi A., Federici L., Formiggini L. 1979, A&A, 35, 
391 

Cid Fernandes R. 1995, PhD Thesis, University of Cambridge 
(available in the WWW: http:/ /www. if.ufrgs.br/~cid) 

Cid Fernandes R., Terlevich R., Aretxaga I. 1996, MNRAS ac- 
cepted 

Cristiani S., Vio R., Andreani P. 1990, AJ, 100, 56 

Cristiani S., Trentini S., La Franca F., Aretxaga I., Andreani P., 

Vio R., Gemmo A. 1996, A&A, 306, 395. 
Di Clemente A., Giallongo E., Natali G., Trevese D., Vagnetti F. 

1996, ApJ, 463, 466 
Edelson R., Krolik J., Pike G. 1990, ApJ, 359, 86 
Francis P. J. 1993, ApJ, 407, 519 

Giallongo E., Trevese D., Vagnetti F. 1991, ApJ, 377, 345 
Hawkins M. R. S. 1996, MNRAS 278, 787. 
Hutchings J. B. 1995, AJ, 110, 994 



